Cotunneling, current blockade, and backaction forces in nanobeams close 

to the Euler buckling instability 
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Single-electron transistors embedded in a vibrating nanoresonator such as a doubly-clamped car- 
bon nanotube exhibit effects stemming from the coupling between electronic and vibrational degrees 
of freedom. In particular, a capacitive electromechanical coupling induces a blockade of the current 
at low bias voltage. It has been recently shown theoretically within a sequential-tunneling approx- 
imation that this current blockade can be enhanced by orders of magnitude when the suspended 
structure is brought to the Euler buckling instability. Here, we investigate the role of cotunneling 
on the predicted enhancement and show that the latter is not suppressed by cotunneling effects. We 
further demonstrate that despite the fact that the current blockade is difficult to measure far from 
the Euler instability, the backaction of the current flow on the nanobeam frequency may be easier 
to observe. 

PACS numbers: 73.63.-b, 85.85. -l-j, 63.22. Gh 



I. INTRODUCTION 

Nanoelectromechanical systems^"'' exhibit effects orig- 
inating from the coupling between electronic and vibra- 
tional degrees of freedom that go beyond the ones en- 
countered in more conventional nanostructures such as 
quantum dots. A prominent example of these new ef- 
fects is the low-bias current blockade induced by a capac- 
itive electromechanical coupling when a single-electron 
transistor is embedded in a vibrating structure.^ This 
phenomenon is the classical counterpart of the Franck- 
Condon blockade in molecular devices^'^ that has been 
observed in suspended carbon nanotubes for high-energy 
vibrational modes. ^ For low-energy (classical) vibrations, 
the blockade is difficult to observe due to the relatively 
weak electromechanical coupling typically encountered in 
experiments. However, it has been reported that this ca- 
pacitive coupling is responsible for the modification of the 
flexural vibration frequency when electrons tunnel on a 
suspended carbon nanotube.^' ^"^ 

In Ref. 11, a way of enhancing the classical current 
blockade by orders of magnitude by exploiting the well- 
known Euler buckling instability^^ has been recently 
proposed: The energy scale associated with the bias 
voltage below which transport is blocked is given by'"' 
E-E = F^/muj'^, where F^ is the electromechanical cou- 
pling resulting from the gate capacitance of the device 
and uj and m the vibrational frequency and mass of its 
mechanical part, respectively. When a lateral compres- 
sive force F is exerted on a nanobeam forming a quantum 
dot coupled through tunnel barriers to source and drain 
reservoirs and capacitively coupled to a gate electrode as 
sketched in Fig. 1, its vibrational frequency uj goes down 
to zero when F reaches the critical force Fc at which 
buckling occurs. ^^ Thus, one expects a large enhance- 
ment of the bias window below which current through 
the device is suppressed. The argument above would in 
principle suggest that the energy scale E^ diverges at the 
instability. However, anharmonicities in the vibrations of 



■^ "7 nanobeam 




FIG. 1. (Color online) Sketch of the setup: a suspended 
doubly-clamped nanobeam (e.g., a carbon nanotube) with de- 
flection amplitude X forming a quantum dot connected via 
tunnel barriers to source and drain electron reservoirs held at 
chemical potentials /il and /ir by the bias voltage V, respec- 
tively. The nanobeam is capacitively coupled to a metallic 
electrode kept at a gate voltage V^, which induces a stochas- 
tic force Fc that attracts the beam towards the gate elec- 
trode whenever the quantum dot is charged. The additional 
compression force F acts on the nanobeam and induces the 
buckling instability. 



the nanobeam, which become crucial close to the Euler 
instability, cut off the apparent divergence. ^^ 

First studied by Leonhard Euler in 1744 while inves- 
tigating the maximal load that a column can sustain,^'' 
the buckling instability is a paradigm for a continuous 
mechanical instability.^^ It has been the subject of var- 
ious experimental studies in micro- and nanomechanical 
systems over the past years. ^'*"^® In particular, buckling 
instabilities have been observed in carbon nanotubes ^'''^^ 
and in Si02 nanobeams^^'^^ where it was demonstrated 
that, despite the small size of these objects, the mechan- 
ical instability can be explained in terms of continuous 
classical elasticity theory. Micromechanical beams be- 
yond the Euler instability were also proposed to serve 
as mechanical memories in Ref. 18. The quantum prop- 
erties of nanobeams close to the Euler instability were 
theoretically studied in Refs. 19-23. In particular, these 



systems were proposed to investigate zero-point fluctua- 
tions of a mechanical mode'^^ or to serve as a mechanical 
qubit.^^'^^ The Euler instability has been recently theo- 
retically considered in nanoeZecfromechanical systems^* 
in the case where the electromechanical coupling F^ is 
negligible with respect to the intrinsic electron-phonon 
coupling. ^^ It was shown in Ref. 24 that this intrinsic 
coupling changes qualitatively the nature of the instabil- 
ity, turning it into a discontinuous one.^^ 

The results of Ref. 11 where the enhancement of the 
current blockade at the Euler instability was predicted 
were obtained in the sequential-tunneling regime of weak 
coupling to the electron reservoirs, i.e., HT <^ k^T with T 
the typical tunneling-induced width (or hybridization) of 
the dot energy levels and T the temperature, where co- 
tunneling currents are exponentially suppressed. In this 
paper, we investigate the role played by cotunneling that 
dominates the current inside the Coulomb diamond and 
which becomes relevant in the opposite regime of strong 
coupling to the leads, i.e., hT ^ k^T. We show that the 
current blockade survives for hybridization to the leads 
smaller than the energy scale E^ set by the bias volt- 
age below which current is suppressed. Moreover, we 
show that the backaction of the current flow on the me- 
chanical part of the device leads to large current-induced 
frequency shifts of the nanobeam that may be easier to 
observe experimentally than the current blockade when 
the nanobeam is far from the mechanical instability. 

The paper is organized as follows: We start by present- 
ing our model in Sec. II. Section III describes the effec- 
tive statistical description in terms of a current-induced 
Langevin process of the nanobeam that we use to solve 
the model. Results for the transport characteristics of 
the device are given in Sec. IV, while Sec. V presents our 
predictions concerning the backaction that the current 
flow exerts on the nanoresonator. We finally conclude in 
Sec. VI. 



II. MODEL 

The model we consider is sketched in Fig. 1: A sus- 
pended nanobeam forms a quantum dot connected to 
source and drain leads via tunnel barriers. A lateral com- 
pression force F is exerted on the nanotube and brings 
it to the Euler buckling instability when F reaches the 
critical force Fc-^^ The total Hamiltonian of the system 
reads 



H = H, 



vib 



H. 



SET 



Hr. 



(1) 



where -ffvib describes the flexural vibrations of the 
nanobeam, Hq^t accounts for the electronic degrees of 
freedom of the single-electron transistor (SET) and H^ 
represents the coupling between the mechanical motion 
of the beam and the electrons residing on the quantum 
dot. 

At sufficiently low temperature, only the lowest-energy 
bending mode of amplitude X is significantly populated. 



and the resulting vibrational Hamiltonian in Eq. (1) takes 
the Landau-Ginzburg form^^'^^"^^ 



Hv 



2rn 



-x^+^x^ 



(2) 



with P the conjugated momentum to X and m the mass 
of the mode. In Eq. (2), 



UJ = wo\/l - 



F 



(3) 



denotes the compression force-dependent vibrational fre- 
quency of the fundamental bending mode of the beam, 
with ujQ the corresponding frequency for vanishing com- 
pression force F. The anharmonicity parameter a > 
stabilizes the mode for F > F^ {uP' < 0) where the 
nanobeam buckles in one of the two metastable posi- 
tions at a: = ±^-TOw2/a. For F < F^ {uj^ > 0), the 
nanobeam stays flat at A = 0. For clamped bound- 
ary conditions, one can express^^ F^ = k{2tt/L)'^, cjq = 
vWf^(27r/L)2, ^ ^ 3crL/8 and a = FcL(7r/2L)'* in 
terms of the bending rigidity k, the length L and the 
linear mass density a of the beam. 

The single-electron transistor is modelled by the 
Hamiltonian 



Hi 



SET 



H, 



dot 



H, 



leads 



Ht, 



with 



-ffdot = (cd - eVg)nd 



(4) 



(5) 



the Hamiltonian describing the single-level quantum dot 
of energy e^, where Vg = CgVg/Cs is the (effective) gate 
voltage which is expressed in terms of the applied gate 
voltage Vg, the capacitance of the gate electrode Cg and 
the total capacitance Ce of the device. The occupation 
operator reads Ud — d^d where d^ (d) creates (annihi- 
lates) an electron on the dot. We consider for simplicity 
spinless electrons, as the inclusion of spin as well as on- 
site Coulomb repulsion in Eq. (5) should not qualitatively 
change our results. The electrons in the left (L) and right 
(R) reservoirs are considered to constitute Fermi liquids. 
The corresponding Hamiltonian reads 



^loads = Xl(^fe ~ ^'')4aCfca' 



(6) 



ka 



where c^^ {cka) creates (annihilates) an electron of mo- 
mentum k with energy e^ in lead a = L,R. The lead a 
is maintained at the chemical potential fia by the bias 
voltage V, such that /^l — A^r = eV- In what follows, we 
assume the voltage drop across the junction to be sym- 
metric, i.e., /XL = — A*R = eV/2. Tunneling between the 
quantum dot and the left and right leads is accounted for 
by the Hamiltonian 



^tun = y^^^'^'^i-a^ + t*ad,'^(^ka), 



ka 



(7) 



with ta the tunnehng amphtude between the dot and lead 
a. 

The flexural vibrations of the nanobeam couple to elec- 
tronic transport due to the gate electrode that exerts the 
electrostatic force —Fc when the quantum dot is charged 
with an extra electron^'^"'^^ (see Fig. 1). The resulting 
coupling Hamiltonian, linear in the deflection X, reads 



H, = F.Xud. 



(8) 



Since the occupation of the quantum dot n^i (which varies 
in our model between and 1) follows a Poissonian statis- 
tics, the coupling (8) results in a stochastic force which is 
exerted on the nanotube whenever the latter is charged 
with an extra electron. Such a coupling induces a block- 
ade of the current at low bias voltage^ which can be dra- 
matically enhanced at the Euler buckling instability.^^ 
Due to the backaction that the current flow exerts on 
the nanobeam, the coupling (8) further induces a renor- 
malization of its vibrational frequency, as recently ob- 
served in suspended carbon nanotubes.^'^^ For quantized 
high-energy vibrational modes, a linear electron-vibron 
coupling similar to Eq. (8), but not induced by the gate 
capacitance, leads to the Franck-Condon blockade. ^"^ 

On top of the electrostatic potential (8) , a second kind 
of coupling between the mechanical and electronic de- 
grees of freedom stems from the electron-phonon coupling 
intrinsic to the tube and is quadratic in X.^^ As shown in 
Ref. 24, such a coupling leads to a current-induced dis- 
continuous Euler instability and to a novel mechanism 
of current blockade at low bias voltage, termed "tricrit- 
ical current blockade". The relative importance of the 
two types of coupling mentioned above is controlled by 
the distance h between the nanobeam and the gate elec- 
trode. Indeed, the intrinsic coupling does not depend on 
h, while the strength of the gate-induced coupling (8) 
depends logarithmically on h.^'^ In the remainder of this 
paper, we assume that the distance h is sufficiently small 
such that the coupling (8) dominates and thus neglect 
the intrinsic electron-phonon coupling. 



III. CURRENT-INDUCED LANGEVIN 
DYNAMICS 

The model encapsulated in the Hamiltonian (1) can 
be solved in an essentially exact way by noticing that 
close to the Euler instabihty, where w — > [see Eq. (3)], 
(i) the vibrational mode is classical, huj <^ k^T with T 
the temperature, and (ii) it is extremely slow as com- 
pared to the electronic dynamics, w ^ F, where F is the 
typical tunneling-induced width (or hybridization) of the 
energy level of the dot. Within a non-equilibrium Born- 
Oppenheimer approximation^^^^^ which exploits the sep- 
aration of timescales between fast electronic and slow vi- 
brational degrees of freedom, and which is asymptotically 
exact at the buckling instability where lu — >■ 0,^^'^"' the 
dynamics of the nanobeam is described by a Langevin 



process with Gaussian white noise, equivalent'^'^ to the 
Fokker-Planck equation 



dtV = - -dxV - F,s{x)dpr 

m 



V{X) 



dp{PV) 



^^l-p 



(9) 



for the probability V{X, P, t) to find the resonator at 
phase-space point {X, P) at time t. The effective force 



Fcff (X) = -muj^X - aX^ - F^no{X) 



(10) 



can be easily inferred from averaging the Heisenberg 
equations of motion for X and P on a timescale long 
as compared to F^^, but short relative to lo^^. In Eq. 
(10), the term proportional to the average dot occu- 
pation for fixed X, no(Ar) = {nd)x, corresponds to a 
conservative force that modifies the dynamical proper- 
ties of the nanobeam in presence of a current flow. In 
Eq. (9), i]{X) and D{X) represent a current-induced 
friction and fluctuation coefficient, respectively. These 
quantities can be expressed through the non-equilibrium 
lesser and greater Green's functions of the dot at flxed 
X, G<{t,X) = i{dU{t))x^ and G>{t,X) = -i{d{t)d^}x, 
respectively. One flnds^^"'^^ 



no{X) = -iJ ^-^G<ie,X), (11a) 

rj{X) = F^ I ^G<{e,X)d,G>{e,X), (lib) 

D{X) = F^J^G<ie,X)G>ie,X), (lie) 



with 



G<(e,X)=y'die-*/''G<(t,X) 



= ih' 



G>{e,X)= Jdtc^^'/^G>it,X) 



(12a) 



= -ih' 



Ea^ail-faie)] 



[e-£4x)]^ + ihr/2r 

In Eq. (12), we defined the Fermi function in lead a 
faie) ^ 



(12b) 



(13) 



exp{[e~ ^ia]/kBT) + 1' 

the effective dot energy for flxed X, 

£d{X) = ed-eVg + F,X, (14) 

and the hybridization F = X^a ^o with F^ = 2T:\ta\'^v /h. 
Here, v is the density of states in the leads, assumed to 
be constant. 

The Fokker-Planck equation (9) can be conveniently 
rewritten using the relevant time, length and energy 



scales of the problem, i.e., the vibrational frequency for 
vanishing compression force coq [see Eq. (3)], the pola- 
ronic shift i — Fc/itiuJq and the corresponding energy 
E^ = Fc£, respectively. Introducing the dimensionless 
position, momentum and time, x — X/£, p — P/mujat^ 
and r = wpi, respectively, one obtains for Eq. (9) 



drV = -pd^r - f,s{x)dpV + i{x)dp{pv) + '^dlv. 

(15) 
The scaled effective force 



/off (a;) = 5x 



ax^ 



no(,xj 



(16) 



is expressed in terms of the reduced compression force 
5 = F/Fc — 1 and the scaled anharmonicity a = a£^/i?£. 
In the low-temperature limit (or equivalently the strong 
coupling limit) k^T <C hV, i.e., the cotunneling regime 
where resonant transport is fully coherent, the occupa- 
tion of the dot for fixed deflection x, the (scaled) current- 
induced dissipation and fluctuation read^^''^^ 



no 



w = 5+E 



— arctan ' ^ 



r/2 



d{x) 



2 — TT 

t^pr^ Y^ To 



T(a^) = -^ J2 



(17) 
^, (18) 



ttT 



Y.^{~^- 



fJ-b) 



a^b 



iia+i 



arctan z 



1 + z' 



r/2 
r/2 



(19) 



/: 



dxdpI{x)Vstix,p), 



(20) 



with 



I{x)^ 



er7L7R 



- arctan 



arctan 



Al + I'g - a; 



r/2 



/ /jR + f g - a: 

V f/2 



(21) 



the quasi-stationary current for fixed x. 

Before we proceed, we notice that for the experiment 
of Steele et al.^ on suspended carbon nanotubes, it was 
estimated in Ref. 11 that E^ ~ 5/ieV and a ~ 10"^", 
yielding F ~ 10, wo — 0.1, and a reduced temperature 
T = k-BT/E^ ~ 1. This experiment is thus in the low- 
temperature adiabatic regime which is investigated here 
where wq ^ T <C F, and where cotunneling might play a 



role on the transport characteristics 
experiments by Lassagne et al.^'^ 



El, 



As for the similar 
~ 5 /ieV and a ~ 



10"^" such that we have F ~ 0.5, Qq - 10"^ and T ~ 
10^. The experiment of Ref. 10 is thus clearly in the 
high-temperature regime Qq <^ T ^ T where sequential 
transport dominates over cotunneling and which has been 
investigated in detail in Ref. 11. 



IV. TRANSPORT CHARACTERISTICS OF THE 
DEVICE 

We now turn to the study of the transport character- 
istics of the device of Fig. 1 with the approach detailed 
in Sec. III. We start in Sec. IV A with a mean-field ap- 
proximation that neglects the current-induced dissipa- 
tion (18) and fluctuation (19). We then consider the full 
non-equilibrium dynamics of the device in Sec. IV B be- 
fore we investigate the effects of a finite quality factor of 
the nanobeam in Sec. IV C. 



A. Mean-field approximation 

It is first instructive to consider the limit wq — > to 
understand the transport characteristics of the system. 
Within this limit, the current-induced dissipation (18) 
and fluctuation (19) can be neglected as compared to 
the effective force (16) entering the Fokker-Planck equa- 
tion (15). In that case, any inflnitesimally small extrinsic 
damping mechanism will, at zero temperature, localize 
the system at a deflection x corresponding to the global 
minimum x^ of the effective potential 



respectively. We introduced the notation 7a ~ Fq/F, 
^g = {e% - ed)/El.,^fli^ = -/iR = v/2 with v = eV/E^, Vcs{x) 
(Do — ^wo/£'e, and F = fiF/Sg. Solving for the station- 
ary solution of the Fokker-Planck equation (15), c^T-'^^st = 
0, one can calculate the average current through the de- 
vice according to 



Sx^ 



ax 



X 

2 



E 

a 

f 



Ja 



(/2a + Vg — x) arctan 
In ( iJla + Vg - x)^ + (f/2)2 



P-a+Vg- X 



F/2 



(22) 



associated to the effective force (16). Within this 
mean-fleld approximation, the stationary solution of the 
Fokker-Planck equation (15) is given by Vst{x,p) — S{x — 
Xni)S{p) such that the current (20) reduces to / = I(xin), 
with X(x) the quasi-stationary current of Eq. (21). 

The mechanism of the classical current blockade is 
encapsulated in the effective potential (22) and in the 
quasi-stationary current (21). Indeed, the effective po- 
tential (22) can, depending on the parameters, show up 
to five metastable minima. These positions x can ei- 
ther correspond to an effective dot energy level (14) lo- 
cated within the bias window, £d{X) < e\V\/2, where 
sequential tunneling dominates and where the current is 
maximal ( "conducting" minima) , either to positions such 
that £d{X) > e\V\/2 where cotunneling dominates and 
the current is algebraically suppressed on a scale given 



by the hybridization F ("blocked" minima). For a given 
gate vohage, the system passes for increasing bias vohage 
from a stable position x^ corresponding to the cotunnel- 



ing (blocked) region to a stable position where sequential 
tunneling dominates (conducting region) , thus defining a 
bias voltage below which current is suppressed. 



1. Limit of vanishing hybridization to the leads. 
We start by considering the limit F — > such that the dot occupation (17) simplifies to 



no(x) 



1, X <v, 

7Le(«) + 7Re(-w), vg-J 

0, 



\v\ 
~2 



\v\ \v\ 

\v\ 
2 ' 



(23) 



where Q{z) is the Heaviside step function. In that case, the dynamical equilibrium condition fcs{x) — and d/cff/dx < 
with /cff (x) the effective force (16) determining the metastable positions of the system can be solved analytically. 
We find that the regions where current can fiow are delineated in the v-Vg plane by straight lines v '^ i2ug. For 
V > 0, and to first order in the small parameter a^/^/|(5| (|(5|/q;^/'^) far (close) to the buckling instability, the apex of 
the resulting conducting region defines a gap A^ below which transport is blocked. 



(24) 



1 

~2S' 








-.5 » 51/3, 


1 
AS' 








(5 » 51/3^ 


1 


lU^' 


-i^ViO- 


-7.-^/^) 


, 1^1 « 51/3. 


251/3-yp 



The gate voltage at which one obtains the minimal threshold (24) reads 

r 7L + 1/2 



26 ' 


7L + 1/2 


4(5 


1 



8(51/3 



7R 



l + 7y^(7R-7L) 



2,5 



1 + 7l^^^(7r-7l) 



-,5 » (51/3, 
(5»di/3, 
\6\ « 51/3. 



(25) 



The case w < is easily obtained by swapping L and R 
in Eqs. (24) and (25). 

The regions in the v-Vg plane where current can flow 
through the nanobeam are shown in Fig. 2 in red/dark 
gray. As it is the case in the sequential-tunneling 
regime, 11 the current blockade characterized by the gap 
(24) increases for increasing compression force 6 < be- 
low the mechanical instability [see Figs. 2(a) and 2(b)]. 
The gap is maximal at the Euler buckling instability 
where (5 = [see Fig. 2(c)] before it decreases above 
the instability [S > 0, see Figs. 2(d) and 2(e)]. Since 
a <C 1,^^ A^ is orders of magnitude larger than away 
from the mechanical instability. Moreover, one observes 
in Fig. 2 a shift of the tip of the conducting regions for 
positive and negative bias voltages which is due to the 



asymmetry of the coupling to the left and right leads 
(7l 7^ 7r)- Indeed, this shift which is maximal at the 
buckling instability [Fig. 2(c)] can be traced back to the 
different occupations of the dot for fixed x [Eq. (23)] for 
positive and negative biases when x is such that the effec- 
tive energy level of the dot (14) belongs to the sequential- 
tunneling region. Such an asymmetry of the conducting 
regions for positive and negative bias voltages has also 
been observed experimentally in suspended carbon nan- 
otubes for quantized vibrational modes^ and is thus not 
specific to the classical regime considered here. 

We conclude this section by noticing that the limit 
F — > presented here is quantitatively similar to the 
zero-temperature mean-field limit considered in Ref. 11 
where sequential transport and symmetric coupling to 
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FIG. 2. (Color online) Conducting regions at mean-field level (in red/dark gray) in the u-Ug plane for F — >■ 0. The red/dark 
gray and white areas correspond to a current |7| = er7L7R and 1 = 0, respectively. The scaled compression force 5 increases 
from (a) to (e). In the figure, 7l = 0.1, 7r = 1 — 7l = 0.9 and a = 10~^. The dashed lines indicate the delimitation of the 
conducting regions for vanishing electromechanical coupling (F^ = 0). 



the leads were assumed, as the effective potential fcff(a;) 
has the same functional form in both cases [compare Eqs. 
(24) and (25) above for 7l = 7r = 1/2 with Eqs. (23) 
and (24) in Ref. 11]. Hence, the enhancement of the 
current blockade is, at mean-field level, not specific to 
the transport regime considered. 



2. Finite hybridization to the leads. 



line), the blockade completely disappears. In that limit, 
one recovers the usual linear-response current in absence 
of electromechanical coupling {F^ = 0), / = 47l7rGoT^, 



with Go 



/27rfi the conductance quantum. The re- 



sults of Fig. 3 show (at mean-field level) that it is ad- 
vantageous to tune the system at the Euler instability, 
where the gap A^ dramatically increases [see Eq. (24)] 
such that the blockade is much more pronounced for a 
given hybridization F to the leads. 



At finite F, we search for the global minimum of the ef- 
fective potential (22) numerically. The resulting average 
current is shown in Fig. 3 as a function of the bias voltage 
for gate voltages corresponding to the apex of the con- 
ducting region, Vg = u™'" [cf. Eq. (25)]. One finds that 
the current behavior, plotted as a function of v/A^, is 
similar for a given r/A„ at the instability [Fig. 3(a)] and 
away from the instability [Fig. 3(b)]. As the hybridiza- 
tion F is increased, the low-bias current blockade becomes 
less pronounced since cotunneling becomes dominant in 
the blocked region {v < Ay). For F ^ A^, (red dotted 



B. Non-equilibrium dynamics 

We now turn to the role played by the current-induced 
dissipation (18) and fiuctuations (19) on the current 
blockade by considering the case Wq t^ 0. To this end, 
we solve for the stationary solution of the Fokker-Planck 
equation (15) numerically as described in Ref. 29. In 
what follows, we focus on the system at the Euler insta- 
bility {5 = 0) where the current blockade is, at mean-field 
level, maximal [cf. Eq. (24)]. Since the stationary solu- 
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FIG. 3. (Color online) Mean- field current 7 as a function of 
bias V scaled by the gap (24) for a gate voltage Vg = u™'" 
corresponding to the apex of the conducting region, see Eq. 
(25). (a) Compression force 5 = corresponding to the Euler 
instability, (b) Compression forces \5\ ^ a^'^ far from the 
buckling instability. In the figure, symmetric coupling to the 
leads is considered (7L = 7r = 1/2). 



tion of the Fokker-Planck equation (15) is, for diQ <C 1, 
independent of the actual value of (Dq (see appendix E in 
Ref. 11 for details), we set wq = 10~^ in what follows. 

Our results for the average current are shown in Fig. 
4(a) for increasing hybridization to the leads T which 
controls the range of the current-induced dissipation (18) 
and fluctuations (19). For F smaller than the energy gap 
A„ (see green squares and orange circles in the figure), 
the current-induced fluctuations have a dramatic effect 
on the I-V characteristics for bias voltages larger than A„ 
[compare with the mean-field results of Fig. 3(a), green 
dashed and orange dashed-dotted lines]. However, for 
F ^ 0.1 At, and v < A„, cotunneling does not suppress 
the current blockade at low bias. We expect that for 
smaller F (which is numerically very difficult to tackle), 
the current blockade should be even more pronounced. 
As close to the mechanical instability, the gap (24) is 
very large, we foresee that the current blockade should 
be clearly visible in an experiment even in the regime of 
strong coupling to the electron reservoirs where cotun- 
neling can be significant. 

The above results are confirmed by the behavior of 
the differential conductance dl/dV shown in Fig. 4(b): 
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FIG. 4. (Color online) (a) Current I including the nonequi- 
librium Langevin dynamics as a function of bias v at the Eu- 
ler buckhng instability {S = 0) for a gate voltage Vg^ = «™'" 
corresponding to the apex of the conducting region, see Eq. 
(25). (b) Corresponding differential conductance d//dV, 
where Go = e^ /2-Kh denotes the conductance quantum. In 
the figure, symmetric coupling to the leads is considered 
(7L = 7R = 1/2), a = 10"'' and Cjq = 10"^ 



For a hybridization to the leads small as compared to 
the energy gap (green squares in the figure) , the differen- 
tial conductance is clearly suppressed at low bias voltage, 
while for larger F, the conductance acquires a finite value 
which is almost constant as a function of the bias voltage 
[orange circles and red triangles in Fig. 4(b)]. 



C. Role of extrinsic dissipation 

Nanoelectromechanical systems are subject to various 
sources of extrinsic dissipation that leads to a finite qual- 
ity factor Q of the nanobeam that we have ignored so far. 
These mechanisms of extrinsic dissipation come from lo- 
calized defects at the surface of the sample, '^'^"'^^ clamp- 
ing and thermoelastic losses,'^'' Ohmic losses due to the 
gate electrode, ^^ etc. Within the Caldeira-Leggett model 
and assuming Ohmic (memory-free) dissipation,"^^ one 
can easily incorporate these mechanisms in the Fokker- 
Planck equation (15) by replacing 7(2;) — > 7(2:) -I- 7c and 
d(x) — > d{x) + 27cT, with 7c = Q^^ the extrinsic damp- 
ing constant (or inverse quality factor) and T — k-eT/E^ 
the (reduced) temperature. 

Numerical results for various values of the ratio 7c/wo 
are shown for the current in Fig. 5(a) and for the cor- 
responding differential conductance in Fig. 5(b). The 
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FIG. 5. (Color online) (a) Average current and (b) differ- 
ential conductance as a function of the bias voltage at the 
Euler buckling instability {5 — 0) for various values of the 
extrinsic damping constant 70. (c) Corresponding effective 
temperature Tiff defined in Eq. (26). In the figure, a = 10~®, 
cDo = 10-^ f/A„ = 0.1, f/A^ = 0.05, 7L = 7R = 1/2, and 
the gate voltage corresponds to the apex of the conducting 
region, see Eq. (25). 



parameter 7e/wo measures the relative importance of 
the two sources of fluctuations, i.e., thermal vs. current- 
induced [cf. Eqs. (18) and (19)], from vanishing extrinsic 
dissipation [green squares in Figs. 5(a) and 5(b)] to purely 
thermal fluctuations (solid hne in the figure), where the 
probability Vst{x,p) entering Eq. (20) obeys a Boltzmann 
distribution at temperature T. We can conclude from the 
results of Figs. 5(a) and 5(b) that the larger the extrinsic 
dissipation (or the smaller the quality factor), the more 
pronounced is the current blockade. As it is the case in 
the sequential-tunneling regime, ^^ the extrinsic dissipa- 
tion helps localizing the system at phase-space points cor- 
responding to the metastable minima of the effective po- 
tential (22), rendering the current blockade sharper. The 
same conclusion holds for quantized vibrations, where 



the Franck-Condon blockade is more pronounced for fast 
equilibration of the vibronic mode.^'^ 

The results for the average current in Fig. 5(a) and 
for the differential conductance in Fig. 5(b) can be inter- 
preted in terms of an effective temperature 



Tcs — 



(rf)/2 + 7er 

(7) + 7c 



(26) 



defined in analogy with the fiuctuation-dissipation 
theorem, '^"^ where 



(7) = / dxdp'^{x)Vst{x,p) 



and 



(d) 



dxdp d{x)Vst{.x,p) 



(27) 



(28) 



are the averages over phase-space of the current-induced 
dissipation (18) and fiuctuation (19), respectively. As 
one can see in Fig. 5(c), T^s almost equals the electronic 
temperature T for inverse quality factors 7c of the or- 
der of ujQ, thus explaining why the current and conduc- 
tance behavior are similar to the case where only thermal 
fluctuations are present [compare red diamonds and ma- 
genta circles with the solid line in Figs. 5(a) and 5(b)]. 
For smaller 70 (larger quality factors Q), the effective 
temperature becomes significantly larger than T and in- 
creases with increasing bias voltage. Hence, the mechan- 
ical system fluctuates more in phase-space and switches 
back and forth between minima of the effective poten- 
tial (22) corresponding to regions where the current at 
fixed x, Eq. (21), is suppressed (cotunneling/ "blocked" 
regions) or enhanced (sequential tunneling/ "conducting" 
regions), thus reducing the average current / and render- 
ing the current blockade less pronounced. 

We conclude this section on the transport character- 
istics of the device by noticing that the features of the 
classical current blockade in the vicinity of a mechanical 
instability are qualitatively the same in the sequential- 
tunneling and cotunneling transport regimes. In the 
former case, the electronic temperature defines the rel- 
evant energy scale below which the current blockade is 
observable, ^^ while in the latter case, it is the hybridiza- 
tion to the electron reservoirs that plays a similar role. 
The enhancement of the current blockade at the Euler 
instability is thus a universal phenomenon that does not 
depend on the transport regime one considers. 



V. BACKACTION OF THE CURRENT FLOW 
ON THE NANORESONATOR 

A characteristics of nanoelectromechanical systems is 
the backaction that the current flow exerts on the me- 
chanical part of the device. The current-induced force 
— no(a;) in the effective force (16) can indeed renormal- 
ize the nanobeam frequency. Since this current-induced 
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FIG. 6. (Color online) (a) Effective frequency ojoir, (b) average occupation of the dot and (c) average current 7 as a function 
of bias voltage v and gate voltage Wg for compression forces far below the Euler instability {—5 ^ a^^^). In the figure, F = A„ 

and 7L = 7R = 1/2. 



force depends on the bias and gate voltages through the 
average occupation of the dot for fixed position x, one 
expects a strong signature of electronic transport on the 
nanobeam frequency. 

In what follows, we quantify the current-induced ef- 
fective frequency Wcff of the nanobeam by the curvature 
of the effective potential (22) at its global minimum a;m, 
'jJo\/v'lii{xm)- Using Eqs. (16) and (17), we have 



Weff 



t^cff == luq. -6 + iaxl 



dnp 



(29) 



Since dno/dvg > 0, this expression shows that the 
nanobeam frequency will be significantly lowered when 
the average number of electrons on the dot varies as a 
function of gate voltage. This is illustrated in Figs. 6(a) 
and 6(b) where we show Woff and the average number of 
electrons on the dot riQ as a function of bias v and gate 
voltage Vg, respectively, for intermediate hybridization to 
the leads F = Ay. For bias voltages v < A^, the aver- 
age occupation of the dot in Fig. 6(b) jumps from to 1 
around a gate voltage corresponding to Eq. (25) and on 
a scale given by the coupling to the leads F, resulting in 
a renormalization of the nanobeam frequency along the 



line Vs 



see Fig. 6(a). For v > Ay, the frequency 



is mostly renormalized for gate voltages corresponding 
to the borders between the sequential tunneling region 
(where uq ~ 1/2) and the cotunneling regions (where 
no — or 1). This is further illustrated in Fig. 6(c) 
which shows the corresponding average current flowing 
through the nanobeam. The results of Fig. 6 are obtained 
for compression forces — 15 3> a^'^, i.e., far below the Eu- 
ler instability. Far above the instability {S ^ <5^^'^), the 
results of Fig. 6 remain unchanged, except that one has 
to replace ujcff/uj by ojeff/v 2|w| in Fig. 6(a) since in the 
buckled state, the vibrational frequency squared in ab- 
sence of a current flow is twice as much as in the flat 
state, see Eq. (29). 

In Fig. 7, we study the influence of the hybridization 
F on the current-induced frequency shift. We show in 
the figure the effective frequency uJcB as a function of the 



bias voltage for gate voltages corresponding to the apex 
of the conducting region [cf. Eq. (25)], for compression 
forces far below the Euler instability (—5 ^ a^^^), and 
for various values of the hybridization F to the leads. It 
is remarkable that for F — > (thick black dotted line in 
Fig. 7), Weft = ^, i-e., there is no effect of the current flow 
on the vibrational frequency of the beam [since in that 
limit, the curvatures of the effective potential (22) is the 
same in all three metastable minima], while the current 
blockade is the most pronounced [see blue solid line in 
Fig. 3(b)]. For F :^ A^,, the effective frequency ujcs/^ — 
1 — 2Ay/TTT is slightly reduced as compared to its value 
without a current flow and is independent of the bias 
voltage (see thick dashed line in Fig. 7) . Only for F ~ A„ 
(see thin lines in Fig. 7) the beam experiences large bias- 
dependent frequency shifts due to the current flow, that 
can reach up to almost 100% (see red solid line in the 
figure) . This is in stark contrast with the behavior of the 
current [Fig. 3(b)], where the low-bias current blockade 
already almost disappears for intermediate values of the 




FIG. 7. (Color online) Effective frequency cjoff of the 
nanobeam far below the Euler instability (—5 ^ a^'"^) as a 
function of bias voltage at the apex of the conducting region, 
Vg — -y^'". In the figure, 7l = 7r = 1/2. 
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hybridization [see orange dashed-dotted line in Fig. 3(b)]. 
The results above show that, although the current 
blockade might not be experimentally accessible far be- 
low (or above) the Euler instability for intermediate hy- 
bridization to the leads, the current-induced frequency 
shift is considerably easier to detect. This was re- 
cently demonstrated experimentally on suspended car- 
bon nanotubes,^'^'^ where despite the fact that the cur- 
rent blockade could not be detected due to the relatively 
weak electromechanical coupling encountered in experi- 
ments, a clear current-induced frequency shift was mea- 
sured. 



VI. CONCLUSION 

In conclusion, we have analyzed the role played by 
cotunneling, relevant in the transport regime of strong 
hybridization to the electron reservoirs, on the low-bias 
current blockade in suspended single-electron transistors. 
We have shown that the dramatic enhancement at the 
Euler buckling instability of the bias voltage below which 
transport is suppressed predicted in Ref. 11 is not ruled 
out by cotunneling currents. The same conclusion holds 
in molecular devices for the Franck-Condon blockade in 



the quantum regime.'^ While in the sequential-tunneling 
regime, the temperature sets the limit of observability of 
the current blockade, it is the hybridization to the leads 
that plays a role qualitatively similar to the tempera- 
ture when transport is fully coherent. The mechanism of 
the classical current blockade and of its enhancement in 
the vicinity of a mechanical instability is thus a univer- 
sal feature of capacitively-coupled nanoelectromechani- 
cal systems, which does not qualitatively depend on the 
transport regime. 

We have further studied the effect of backaction forces 
on the mechanical properties of the nanobeam. We have 
shown that the mechanical frequency can be strongly 
renormalized as a function of the bias voltage due to the 
force that is exerted by the current flow on the nanobeam. 
This effect might be considerably easier to observe than 
the current blockade when the system is far from the 
Euler buckling instability. 
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